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Abstract. The modified Bernoulli numbers 

V 2r y n + r 

7' — 

introduced by D. Zagier in 1998 were recently extended to the polynomial case 
by replacing Br by the Bernoulli polynomials Br{x). Arithmetic properties of 
the coefficients of these polynomials are established. In particular, the 2-adic 
valuation of the modified Bernoulli numbers is determined. A variety of ana- 
lytic, umbral, and asymptotic methods is used to analyze these polynomials. 



1. Introduction 
The Bernoulli numbers J5„, defined by the generating function 

(1-1) ^ = E^A' 

n—O 

were extended by D. Zagier |17) with the introduction of the so-called modified 
Bernoulli numbers i?* defined by 

(1.2) b: = ±^-+^^ 



r=0 



2r J n + r 



Note that B'q is undefined. Arithmetic properties of B2n {Bi = — ^ and B2n+i = 0, 
for n > 0), include the von Staudt-Clausen theorem which states that, for n > 0, 



(1.3) B^n^- J2 - 



(p-l)|2„ 
p prime 



mod 1. 

P 



It follows that the denominator of B2n is the product of all primes p such that 
p — I divides 2?!. On the other hand, the numerators of B2n are still a mysterious 
sequence. 

The definition (|1.2p shows that B* is a rational number. Write it in reduced 
form and define 

(1.4) an = denom(B*). 
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Zagier JJj showed that 

(1.5) 

satisfies 

(1.6) E 



Bn = 2nB^ — Bn 



E - 



mod 1, {n > 0, n even) 



(p+l)|r. 
p prime 



that imphes 
(1.7) 



2nBl 



(p-i)l" 

p prime 



P 



(P+l)l 
p primi 



1 



mod 1, (n > 0). 



This statement shows that if p is a prime dividing a„ (defined in (|1.4p ). then at 
least one of p, p — 1 and p + 1 divides n. In particular, all prime factors p of a„ 
satisfy p < n + 1. In fact, from computing the first 1000 terms, it appears that, 
conjecturally, the following stronger statement is true: if p is a prime dividing q;„, 
then p + 1 or p — 1 divides n. 

The first few values of the sequence {B*} are 

3 1 1 27 1 29 1 451 1 65 

4' 24' ^4' ^80 '^4' ^1260' 4' 1120' 4' ~ 264' ' ' ' 

Our particular interest will be in the 2-adic properties of this sequence and the 
2-adic valuation of B* will be worked out completely. A guiding question motivated 
by the first few terms as above is: 

Question 1.1. Is the denominator a„ always divisible hy A? 

This basic question will become particularly relevant when considering the cor- 
responding modifications of Bernoulli polynomials. This is addressed at the end of 
this introduction. 

It turns out that a2n+i = 4, so only even indices need to be considered. The 
first few values of \a2n are given by 

(1.8) 6, 20, 315, 280, 66, 3003, 78, 9520, 305235, 20900, 138, 19734, 6, 7540, . . . 

This sequence has been recently added to OEIS (the database created by N. 
Sloane) as entry A216912. The next figure shows the 2-adic valuation of a2n] that 
is, the highest power of 2 that divides a2n- 





60 80 



Figure 1 . Power of 2 that divides denominator of 
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Symbolic computations lead us to discover the next result. In particular, this 
answers Question ll.ll in the afBrmative. 



Theorem 1.2. For n > 0, 

{1 if n = 6 mod 12, 
2 if n = mod 12, 
otherwise. 

In particular, _B*, the denominator of an, is divisible by 4. 

Note that this may be rephrased in the following way: The 2-adic valuations 
V2{SnB2n) form a periodic sequence of period 6 with values 

(1.9) {0,0,1,0,0,2}. 

This is an unexpected variation on the period 6 theme: D. Zagier proved that the 
sequence {B2n+i\ is 6-periodic. 

The modified Bernoulli numbers -B* were extended in [B] to the Zagier polyno- 
mials defined by 



(1.10) Bl{x) = Y. 



r=0 



n + r\ Br[x) 



2r J n + r 



so that = -B*(0). The first few are: 

i(2a; + 3), ^ {6x^ + 18x + l) , ^{2x + 3) (x^ + 3x - l) , 

— (lOx'^ + 60x^ + 90x^ - 27) , — (2a; + 3) (3x'^ + 18x^ + 23x'^ - 12a; - 5) , . . . 
80 60 

In analogy to a„ in (|1.4p . define, for j E Z, 

(1.11) a„,j = denom(B*(j)). 

It is shown in Lemma 13.21 of Section [3] that, under the assumption that 4 divides 
a„, the denominators a„ equal q;„ for any j g Z. Combining this with Theorem 
11.21 one obtains: 

Theorem 1.3. The denominator a^.j = denom(i3* (j)) does not depend on the 
value j G Z. 

Special values of B^[x) present interesting arithmetic properties. The relation 

(1.12) Bl{x + 1) = Bl{x) + ic/„_i (I + l) , 

relating iJ* to the Chebyshev polynomial of the second kind, appears as Lemma 
10.2 in "6^. In particular, this shows the identity 

n 

(1-13) Bl{l)=Bl + -. 

On the other hand, the values i?*(— 1) are connected to the asymptotic expansion 
of the function 

(1.14) 1/(2) =logz + V U+i 
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at z — >■ 0. Here, ip{z) is the digamma function 

(1-15) m - 

the logarithmic derivative of the gamma function. The proof of the next statement 
appears in Section [71 

Theorem 1.4. Define the numbers Vn by the asymptotic expansion 

oo 

(1.16) y(z)^^i;„z". 

Then w„ = -2B;(-1). 

The value W2n-i ~ (^l)"/2 is simple to obtain, but 



(1.17) V2n = (-1) 



n+1 



n ' [ n-k I 2k 

k=l ^ ^ 



requires further work. 

A second motivation for considering the sequence comes from the natural 
interest in the sequence {B\^. The established fact that {i32ji+i} is 6-periodic has 
no obvious analog for the even indices. It turns out that the function Viz^ satisfies 



oo -. 

(1.18) E^2V'" = -9^W-^ 



2 ' ' 4 

n=l 



1 2(1 



1 1 



thus connecting i?2„ and w„. 

A variety of expressions for the coefficients Vn are provided. Section [6] gives 
one using the umbral method and Section [7] exploits a relation between the Zagier 
polynomials _B* and the Chebyshev polynomials Un{x) to determine Vn- A direct 
asymptotic method is used in Section |5] and Section [S] presents a family of poly- 
nomials that determine w„. The classical integral representation of the digamma 
function is used in Section [TUl the formula of Faa di Bruno to differentiate compo- 
sitions is used in Section [TT] and, finally, a recurrence for Vn is analyzed in Section 
[TUby the WZ-method [H]. 

2. The 2-adic valuation of B* 



The goal of this section is to establish Theorem 11.21 which determines the 2-adic 
valuation of the sequence 

The strategy employed here is as follows. It is a consequence of the von Staudt- 
Clausen congruence that the Bernoulli numbers 2i?„ are 2-integral. From this one 
may conclude that the rational numbers 47T.J5* are 2-integral as well. In particu- 
lar, these numbers can be reduced modulo powers of 2 to determine their 2-adic 
valuation. Here, it will be sufficient to reduce them modulo 8. To begin with, the 
classical Bernoulli numbers are reduced modulo 8. 

Proposition 2.1. The following congruences hold modulo 8: 

1 if k even, 



2Bq = 2, 2i?2 = 3, 2i?2fe 
with k > 1. 



5 if k odd, 
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Proof. The von Staudt-Clausen theorem states that 

(2.1) pB2fc=p-lmod/+i 

for p prime, k > 2 and when (p — l)p^ divides 2k; see [13], formula 24.10.2 on page 
593. Now take p = 2 and £ = 1. Then for > 2 it follows that 2B2k = 1 mod 4. 
Therefore 2B2k = 1 or 5 mod 8. In the case k is even, one may take £ = 2, since 
then (p — l)p^ = 4 divides 2k. Therefore 

(2.2) 2B2k = 1 mod 8. 

A different proof of this fact appears in . The identity established there is 

(2.3) 2B2k = 1 mod 2'+^ 

where 2'' is the highest power of 2 that divides 2k. In particular, for k even, r > 2 
and the result follows. 

The case k odd requires a different approach. 



Let Urn be the numerator and Vm the denominator of Bm, so that Bm = Um/Vm 
and {Um,Vm) = 1, Vm > 0. Voronoi's congruence [TTl Proposition 15.2.3] states 
that, if m > 2 is even and a, n are positive integers with (a, n) = 1, then 



ri-l 



(a" - l)C/„, EE ma"- V,„ ^ j"'-' 



n 



mod ; 



As usual, [x] refers to the greatest integer less than or equal to x. It follows from 
the von Staudt-Clausen congruence that 2B2m has 2-adic valuation for to > 0, 
so that they are 2-integral. Voronoi's congruence with a = 3 and n = 64 therefore 
yields 



63 



(3™-l)2S„ = 2TO3™-i^j'"-i 



64 



mod 64. 



One easily checks that, for even m, 3™ — 1 = 4to modulo 64. Similarly, after 
checking finitely many cases, for m = 2 modulo 4 with to > 6, 



63 



i=i 



64 



= 42 mod 64. 



Combining these, one finds, for to = 2 modulo 4 with m > 6, 



TO 

2Bm 



TO 

5— mod 8. 



Hence, if to = 2fc with fc > 3 odd, then 2i?„, = 5 modulo 8. 



□ 



Further basic ingredients are the following generating functions. 



6 



M. COFFEY FT AL. 



Proposition 2.2. The following generating functions admit rational closed- forms: 
(2.4) 2 + ^a;"^f ) ^TTI " 1 - 3x + ' 



n=l fe=0 
L"/2J 



2 « Y-V" + 2n _ (1 - 2x) (2 - 2x 

r! = l fc=0 ^ ^ 
n=l fc=0 ^ ^ 



4fc )n + 2k (l-.T + a;2)(l-3x + x2)' 
n + 2fc\ 2n 2-6x + 7x^ - 2x^ 



Proof. These readily follow from the generating function for T„(a;), the Chebyshev 
polynomials of the first kind, given by 



(2.5) Y.T„{x)t-^- 
and from the fact 



1- xt 



^ - 2xt + t2 

n=0 



(2.6) y " + n^^iT„f^ + i 

^ ' ^\ 2r 1 n + r n"V2 

proved as Lemma 9.1 in [5]. □ 



Equipped as such, a proof of Theorem 11.21 is given next. The statement of this 
theorem is repeated for the convenience of the reader. 

Theorem 2.3. For n> 0, 

{1 if n = 6 mod 12, 
2 ifn = mod 12, 
otherwise. 

Proof. It is convenient to remark at the beginning that the case of odd n is simple 
and is a consequence of Zagier's result on the periodicity of the sequence ^2™+!- 

Working modulo 8 and using Proposition l2.1[ it follows that 2Bo = 2, 2Bi = — 1, 
2B2 = 3 and for k > 1, 

2B2k = 3-2i-lf . 

Note that ("2'^'') is an integer. Thus it follows from (ITT^ that 4nB* is a 
2-adic integer. For n > 1, these numbers reduce modulo 8 to 

k=0 



k=0 



= -n^ + 2-n 



2k J n + k 

L»/2J 

4fc J n + 2k 



y r+^'']^2B,k 

^ \ Ah T, 



n+l 
3 



- - - + 2k\ 2n 



E [3-2(-l)* 



fe=0 



Ak n + 2k' 
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where in the second equahty, the term conies from the contribution of Bi — 
— 1/2, the only nonzero BernouUi number of odd index. Also, for the final congru- 
ence, adjusting for the k = and fc = 1 cases in which 2Bo = 2 and 2B2 = 1/3 = 3 
respectively, produces the extra terms 



n\ 2n 




(2Bo - 1) 



n + 2 
4 



2n 



{2B2 ~5) = 2-n 



n + 1 
3 



Using Proposition 12.21 modulo 8 now gives 

2 X {I + x) {1 + x'^) 
(1-x)^ 
-6x + 7x'^ - 2x^ 



n=l 



(1 - 2x) {2-2x + x^) 
'(1 - a; + a;2) (1 - 3x + x"^) 



- 2- 



l-Ax + lx'^- 4x3 _|_ ' 
where it is readily verified that the right-hand side is a rational function whose 
coefficients modulo 8 are periodic with period 24. The even part simplifies to 



^8ni?*„x2"^ 



xii 



6x~ 



■ 3x^ + Ax^ 



n=l 



This implies 

{8nB*„) = 

which proves the claim. 



if (n,3) = l, 

if n = 3 mod 6, 

if n = mod 6, 



□ 



The goal of this section is to establish Theorem ll.3l It states that the denomi- 
nator of B^{j) does not depend on j € Z. The proof begins with the identity 



3. The denominators of B*{j) 



(3.1) 



B:{x + l) = B*ix) + ^Un 



appearing as Lemma 10.2 in |6j which establishes a relation between the Zagier 
polynomials and the Chebyshev polynomials of the second kind Un{x). 

Lemma 3.1. For every half-integer x, the numbers Un{x) are integers. 

Proof. This is clear upon using the determinant representation 

2x 1 



(3.2) 



Unix) 



1 2x 



1 

2x 



1 

for the Chebyshev polynomial. To verify p.2p denote the determinant by Dn{x). 
By expansion by minors, it follows that Dn+i{x) — 2xDn{x) — Dn~i{x). The same 
recurrence is satisfied by Un{x) and a direct computation gives Dn{x) = Un{x) for 
n = 1, 2. Thus, Unix) = Dnix) for all n G N. 

An alternative proof employs the generating function of the Unix) polynomials 

(3.3) J2Ukix)t''= ^ 



fc>0 



l-2xt + t'^' 
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Choosing x — ^ with p integer, it follows that 

(3.4) y Urn ( ^) = , ^ = , r V = F t'iP - t)" 

\ ) '^\2J 1-pt + t^ l-tip-t) ^ "-^ ' 

m>0 ^ ' k>0 

since by choosing t small enough, \t{p — t)\ < 1. The coefficient of in this sum, 
which is Um (f), is clearly an integer. □ 

Lemma 3.2. The denominator of B*(j) is independent of j € Z. In other words, 
for all j G Z, 

(3.5) denom B * ( j ) = denom _B * . 

Proof. Assume j > 0. It is a consequence of Theorem II . 2 1 that the denominator of 
S* is divisible by 4, and thus is At for some i € Z. 

Assume, therefore, by induction that the denominator of B^{j) is At as well; that 
is, in reduced form 

(3.6) s:(j) = 

with X ~ x{j) an odd integer. The identity p.ip coupled with Lemma 13.21 gives 
^\ / ■ -, X 2^ ^ 2; + 2wt 

(3.7) ^"(^ + l)-4^+2=^-4^' 

with w e Z. The last fraction in (13. 7p is also in reduced form. Indeed, the numerator 
is odd so there is no cancellation of the factor 4 and if p is an odd prime that divides 
both X + 2wt and At, then it divides gcd{x,t) — 1. Therefore B^{j + 1) also has 
denominator At, the denominator of B*. This proof easily adapts to the case when 
j is negative. □ 

4. An asymptotic expansion related to the numbers -B* 
The generating function 

= -2 log^ - 2'^' {z+l/z-1- x) 
appears as Theorem 5.1 of f6l. Here 'ip{z) is the digamma function 

(4.1) i^{z) = -f^, 

r(z) 

the logarithmic derivative of the gamma function. The asymptotic expansion for 
the auxiliary function 

(4.2) l/(z) =.logz + V U+- 



as z — > in the form 



(4.3) V{z) 



00 



Vr,Z 



will yield a relation between the numbers -B* and the sequence Vn in (14. 3p . 
The value of a2n+i has been established in [5]. 

Theorem 4.1. For j G Z, the coefficients 4-B2n+i(i) ore odd integers. This gives 
(4.4) a2„+i = 4. 
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The generating function for the much more involved case of a2n is 

2 log Z — -il> I ^ ^ -7/j \ y. 



E BLUy = --logz--ij(z + - + 2 + j)~-i; 



z + - - 1 - J 

z 



This was given in CoroUary 5.3 of 6| and can be converted to 



1 



log z — —ip [ z + 



2 ° 2 
4 



1 



1 



r=0 



using 
(4.5) 



rz + 1 z^ — rz + 1 



k-l 



4:{Z^ + 1) 



■"i — ^ n 4- 



r=0 



u + r 



Now use the function V{z) defined in (|4.2I) to obtain 



(4.6) 



n=l 



Z^ + 1 



1 - Z*^ 



This identity shows that Question 11.11 is indeed equivalent to the rational num- 
bers V2n having even denominators. 

A direct symbolic computation gives the values of the first few Vn as 

r 1 11 1 13 1 29 1 109 1 67 1 65711 
^ ' ' t ' ^2' 12' 2' ~40' ~2' 630' 2' 560' "2'^T32' 2' 6006 J ' 

This data suggests that \vn\ — 1/2 for n odd but no simple pattern is observed for 
n even. 

5. The use of bounds on 'ip{z) 

The first approach to the computation of the coefficients w„ is to use bounds for 
the digamma function V-'(^) ^-nd its derivatives that exist in the literature. This 
process succeeds only for small values of n. 



Proposition 5.1. The function V{z) satisfies 
(5.1) 

that is, vq — 0. 
Proof. The inequality 
(5.2) 

was established by H. Alzer [2]. This gives 

1) - < V{z) < log(z2 + 1) 



lim V(z) = 0, 

2^0 + 



< logz - i;{z) < - 
2z z 



(5.3) log(z2 

and the result follo\ 
and B. Guo [TO] to 



Z2 + 1 ' ' • ^ 2(22 + 1) 

and the result follows from here. The inequality (|5.2p has been improved by F. Qi 



(5.4) 



log 



(^z + - i < V(^) < log (z + e "')-^- 
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□ 

The next statement shows the computation of vi . It requires sharper bounds on 
the derivative ip'{x). The proof presented below should be seen as a sign that a 
different procedure is desirable for the evaluation of general w„. 

Proposition 5.2. The function V{z) satisfies 
(5.5) linirW^-^, 
that is, vi — ^1/2. 
Proof. The inequalities 

(5-6) ^ + J^<(-irV«W<^ + ^, for.>0, 
are established in [5]. In the special case k = 1 they produce 

5.7 - + -^ z < - + -. 

It turns out that the lower bound gives a sharp result for V'{z) as z 0+. Indeed, 
(5.8) r(z)= z+- +-< 



V ^ y V zj Z 2(l + z2)2 

The reader should check that the upper bound does not give useful information. 
Instead the inequality 

(5.9) V'(^) < e'/" - 1, 

established in [9], is used to produce 



(5.10) y'(z) > 



1 



exp 1 -r—r ) - 1 



,z^ + 1. 

The result now follows by letting z — > in (|5.8p and (|5.10p . □ 

The computation of w„ by this procedure requires bounds on all derivatives of 
ip{x). The examples discussed above shows that this is not an efficient procedure. 
The next section presents an alternative. 

6. The computation of u„ by umbral calculus 

The goal of this section is to compute the coefficients w„ in the expansion (|4.3|) by 
the techniques of umbral calculus. The reader is referred to [G] for an introduction 
to these techniques and for the statements used in this section. 

Introduce the auxiliary function 

(6.1) i^(x) = V +loga;, 
for a; > and observe that 

(6.2) v{z) = F(^-^^+log{z' + l). 
Theorem 6.1. The function F{x) admits the asymptotic expansion 

(6.3) F{x) ^ \ 

n=l 
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Proof. The integral representation 

(6.4) ^(z) = log z + e-*- Q - Y^-n) dt 
produces 

(6.5) Fi.)^fe-^^^(l-^)dt. 
Set s — t/x to obtain 

The generating function for the Bernouhi polynomials 

(6-7) ^^^E^"(^)- 

e"^ — 1 — ' n! 

yields 



2-^ Til 



s '■ — ' n\ 

Tl=l 



^1 ^0 



The result now follows from i?n(l) = (— l)"i3„. □ 
Note 6.2. The asymptotic behavior 

(6.8) |S2„|~4V^ — 

VTre/ 

shows that the series in (j6.3p does not converge for x 7^ 0. 



The result in Theorem l6.1l is now transformed using the umbral method described 
in [6]. The essential point is the introduction of an umbra *B for the Bernoulli 
polynomials Bn{x) by the generating function 

(6.9) eval{exp(i<8(a;))} = 4^ 

— 1 

The rules eval (*B") = i3„ and eval (*B(a:)) — eval {x + *B} are useful in converting 
identities involving Bernoulli polynomials. 

Theorem 6.3. The coefficients Vn in the expansion (|4.3p are given by 

(6.10) Vn - V i-ir-''+'^^Bn-2k. 

k=0 

Proof. The result of Theorem 16.11 can be written as 

(6.11) F{x) = — '- (x«B)" 

n—l 

= log(l+x<B). 
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Then 

V{x) = F(^-^)+log(x2 + l) 

= eval (log (1 + / ] + \og{x^ + 1) 



eval(log {x^ + 1 + a;*8)) 
eval ^ — x''{x + <B)" 

\r=l ^ 

'- X'-Brix) 



r=l fc=0 ^ ^ 

Now let n = r + fc and invert the order of summation to obtain the result. □ 

Separating the expression for the coefficients w„ given in ()6.10p according to the 
parity of n, simplifies the result. 

Corollary 6.4. The coefficients w„ in (|4.3p are given by 

(-1)" 



(6.12) V2n-1 = 



(6.13) V2n = (-1) 



n+1 



1 " 

n ^ — ' \ n 

k=l 



n + k-l\ B- 



>2k 
'2k 



7. Properties of Zagier polynomials give the expression for w„ 

This section presents a proof of the expressions for w„ given in CoroUarv 16.41 bv 
using properties of the Zagier polynomials established in [B]. 

Theorem 5.1 in [6] gives the generating function of the Zagier polynomials 

(7.1) gB*(,),"^_i2|f + 

that for a; = —1 yields 

n=l ^ 

Comparing with the asyniptotics for V{z) given in (j4.3l) gives the next statement. 
Proposition 7.1. The coefficients u„ are given by 

(7.3) = -2i?:(-l). 

To obtain an expression for i?*(— 1) use p.ip with n replaced by 2n + 1 and 
X = — \. It follows that 

(7.4) i?2*„+l(-l) = 52*„+i(0) - \u2n (\ 
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The reduction of this expression uses Theorem 10.1 in [6 in the form 



U2n 



r-0 ^ ' / ' 

which in the special case x = produces 

f-D" 1 ^l 

(7.6) B2Vi(0)-^+2^2„(^2 

using C/2n(0) = (—1)". Inserting in (17.41) gives the resuh for odd index. 

In the case of even index, the proof starts with the reflection symmetry of the 
Zagier polynomials 

(7.7) b:{-x - 3) = {^irB:{x) 

(given as Theorem 11.1 in [6]) which in the special case x = —2 gives 

(7.8) -2B*„(-l) = -2i?2*„(-2). 

To obtain the expression for V2m use the identity (10.10) in [5] 



2r I n + r 

in the special case x — 0. This gives the values of w„ stated in CoroUarv 16.41 Thus 
([6T3| and dLH]) imply ([731 . 



8. Calculation of Vn by an asymptotic method 

The goal of this section is to derive the formula for by a direct asymptotic 
expansion of the digamma function: 



^ oo 

(8.1) ^^,)^iogz---J2 
Start with 



2z ^ 2fcz2fe 



B2k 

as z — )■ OO. 



5.2) y(^) = ^(£!±l') _iogf£!±i)+log(z2 + l) 



and use (18.11) to obtain 



( z 



^'w ~ -Ef brr) +'°^<- + " 



fe=i 

2fe 



2(z2 + l) 



f- (-lr+^ 2n ^^n,2„+l ^ i?,, ^ (- 1)^ (£ - 1 + 2fc) ! 

~ ^ n 2^^ ' ^ 2k ^ t\{2k-V 

n=l n=0 k=l e=0 ^ 

The coefficient of the odd powers of z can be read immediately. Indeed, 

(-1)" 

(8.3) V2n-1 = ^7^. 
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This is (|6.12p . To obtain the expression for the even powers, observe that 

^B^2k^ (-1)^(^-1 + 2fc)! _ ^f,y ( (^ + k-l\B^\ 2. 
2^ 2k ^ l\{2k - 1)! ~ M M 2fc - 1 ; 2fc I 

This gives 

k=l ^ ^ 

This is equivalent to (|6.13p and also to (I7.9P with a; = 0. 

An expression for V2n in terms of Chebyshev polynomials in given next. 
Proposition 8.1. Let Tn{x) he the Chebyshev polynomial of the first kind. Then 

(8.5) W2n = -em/ ^^T2„ ^y^^ = -^Tn 
Proof. Lemma 9.2 in |^ established the representation 

(8.6) B;(^)^eval(iT„(^±^ 
The relation (|7.8p gives V2n = — 2i?2„(— 2) so that 

(8.7) V2n = -eval (-T2n 



The result now follows from the identity T2n{x) — r„(2a;^ — 1) for Chebyshev 
polynomials; see [3 , 7.210 formula 7 on page 550. □ 

9. The asymptotics of 'iJj{z) and its derivatives 
The coefhcients t;„ in the expansion (14. 3p are now evaluated from the expression 

1 f d \ ^ 

(9.1) Vn = ^ lim — V{z). 

n\ 2^0 \ dz J 

The next theorem shows existence of a sequence of polynomials ylj.„(z) that give 
the desired formula for derivatives of ^(z). Theorem 19.21 presented below provides 
an explicit form of these polynomials. 

Theorem 9.1. Let n G N. Then there are polynomials Aj^n{z), with 1 < j < n 
such that 

(9.2) z'- i^-j V{z) = {-ir-\n ~ l)!z" + A„„(z)^, (z + 1/z). 
The polynomials Aj.„(z) satisfy the recurrences 

^n+l,n+l(^) = (z^ - l)A„,„(z), 

^j.n+i(^) = -2nzAj-„(z) + z2A;.„(z) + (z2-l)Aj-_i^„(z) /or2<j<n, 
Ax^n+\{z) ^ ^2nzAi^n{z) + z^A[Jz), 
and the initial condition 

Ai^iiz) ^ z^ - I. 

The degree of Aj_„(z) is n + j — 2 ifl<j<n— 1 and 2n for j = n. 
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Proof. The term (—1)" ^{n ~ l)lz " arises from the n-th derivative of \ogz. To 
obtain the recurrences, simply observe that 



n+l 



dz 



and compare the coefficients of + l/^)- The statement about the degree of 
Aj^n{z) is obtained directly from the recurrence. □ 

The next theorem gives an explicit form of the polynomials Aj^n{z). The authors 
wish to thank C. Koutschan who used his symbolic package to solve the recurrences 
in Theorem 19. II 

Theorem 9.2. The polynomials Aj.„(z) are given by 
(9.3) A„,„(z) = (z^ - 1)" 

and for 1 < j < n, 



(9.4) 



A,,„(z) = (-l)"^z"-^-^(-l) 

^ ' r=0 



n — \ — r\ f j 
n-j 



Proof. Simply check that the form stated in this theorem satisfies the recurrence 

□ 



given in Theorem 19. II 
Note 9.3. The package of C. Koutschan actually gives the form 



(9.5) A,,„(z) = 



n-j\/n 
J - V \J. 



{n-jy.2Fi {l^j.-j;l-n;z^) 



The hypergeometric representation of the Jacobi polynomials 



(9.6) 



2F1 —m, TO + a + /3+l;a + l; 



1-x 



shows that 
(9.7) 

Note 9.4. The coefficients Vn are now obtained from ()9.ip and the identity 
(9.8) 

This employs the expansion 



] 



(9.9) i^^^\z)^i;,{z)^{-iy-^ 



Z3 



2zi+ 



k=l 



( 2fc + j-l)! 
(2fc)!z2'=+J 



(that appears as 6.4.11 in Tl). The polygamma function, which appear differenti- 
ating (|4.2p to obtain (|9.1I) . has argument z + 1/z. Thus (|9.9p is used in the form 



(j - ly.z^ Jlz3+^ B2k{2k + J - l)!z2^-+^- 

" ■ 1)^+1 (2fc)!(z2 + l)2fc+J 



- ^^'y^' (.2+1). ■ 2(.2+i) 

(z2 + l)j 2^ k\ 



k=0 



(z2 + l)fc 
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as Z 0. 

Proposition 9.5. The asymptotic expansion 



2 

r=0 



holds as z ^ 0. 



i-kr. {k + j + e-iy. 



z^' 



A direct non-illuminating computation of the expansion in (j9.8l) gives the values 
of Vn in Theorem 16.41 Given the fact that other proofs of this result have been 
provided, the long but elementary details are omitted. 

10. Calculation of Vn via integral representations and the Faa di 

Bruno formula 

This section employs the integral representation 

/■°° / I 1 

(10.1) iljix) =loga; + 



t 1 - e 



of the digamma function, given as entry 8.361.8 in f7|, to obtain the values of Vn 
given in Corollarv 16.41 



Lemma 10.1. The Junction V{z) in (|4.2p is expressed as 

(10.2) V{z) - log(z2 + 1) + _ dt. 

The representation (jl0.2l) reduces the computation of w„ to the asymptotic ex- 
pansion of 

(10.3) W{z) = Q - ^ dt. 
Indeed, if 

(10.4) V{z) - E i;„z" and W{z) - ^ w„ 



OO CO 

ri-1 , 



then V2n-i = W2n-i and V2n = + (^1)" /n. The next lemma is preliminary 
to the computation of this expansion. 



\ 2n OO / , ; 1 

z \ , ,1, .„ I n + k — I 



Lemma 10.2 

, -V(-i)''"", 

^z^ + lj ' \ k-n 

Proof. Use the binomial series 



^2k 



Z 



to find 

/ \ 2n OO / OO / ^ 

„2fe 



2—0 fc— n 
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Now use the elementary identity 

-271 \ ^ ,^.k~nfn + k- 1 
k — n J \ k — n 

to obtain the result. □ 

To find the asymptotic expansion of the function W{z) defined in (|10.3p . let 
s = z/ (z^ + 1), and use the change of variable x — t/s to get 



Jo ^t^, ^! 



The infinite series is not uniformly convergent as z — )■ 0, and interchanging the 
sum with the integral does not provide a convergent series. But the resulting series 
(with radius of convergence zero) will be the asymptotic expansion of W{z): 



oo 



W{z) ^ -Y^^^s'^ r :,^^-^e--dx 



n=l 



oo 



n—l n— 1 

2(z2 + l) 2^ 2n U2 + I 

The expression for the coefficients Wn corresponding to (j6.13l) now follows from 
Lemma 110.21 

An alternative approach based on the integral representation 110.11 uses the Faa 
di Bruno formula and the partial Bell polynomials. Write 

1 1 

1 ^ 1 - e- 



i^[x) = i^{x) -\ogx^ I ( - - — ] e--'dt, 



so that W{z) = ■ilj{h{z)) with h{z) = z + l/z and 



e-'^'dt 



Define 



(10.5) hiz) = - ^^-"^'^'^'^dt = V;W(M^)). 

The partial Bell polynomial Bn^k in the n — k + I variables xi, . . . , Xn^k+i is 
defined by 



Bn,k{xi,...,Xn-k+l)= V . , . , r (77) (^) 



.Tiyi 1^X2^^ I Xn-k+l 



. ^ii\j2\---jn-k+i\\VJ \2\) \{n-k+l)\, 

a{n,k) 
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where the sum is over the set a{n, k) of all non-negative integer sequences ji , j2 , • • ■ , jn-k+i 
such that 

ji+j2-\ \-jn-k+i=k and ji + 2j2 H \- {n - k + l)jn-k+i = n. 

The Faa di Bruno formula for the n-th derivative of the composition W = ip o h 
is then expressed as 

n 
n 

(10.6) = ^/fc(z)i?„,fc(/i'(z),... ,/i("-^-+i)(z)). 

The next lemma provides some results on the partial Bell polynomials. A useful 
reference is 5 , page 133-137. 

Lemma 10.3. 

(10.7) BnM{xi,st^X2, st^X3,st'^X4, • • • ) = s''f^Bn,k [^,X2,Xz, ' ' ' j 

(10.8) B„,fe(xi, X2, . . .) = l4Bn^k,k~e (y , y, • ■ •) 

(10.9) S„,.(l!,2-3- 

Proof. The proof of (|10.7p follows easily from the definition, noting that 3j2 -I- 4^3 -I- 
■■■ =n + k-2ji. Formula (fTOSl) is entry [3.1] on [51, and (fTa9| is entry [3.h]. □ 

Lemma 10.4. The partial Bell polynomials satisfy 

(10.10) B„,.(/^'(z),... ,h^-'^'\^))-^-^i:l 

1=0 

Proof Note that h'{z) = 1 - z'^, and h'-'\z) = {-lyUz-'-'^ for i > 1. Hence the 
result easily follows from (with s = -l/z,t ^ l/z), (fTOTS)) and ((TOJ)) . □ 

The next result expresses the integrals Ik{z) defined in (jlO.Sp in terms of the 
Hurwitz zeta function 

CXD 

Proposition 10.5. The integral Ik{z) is given by 

^'^'^ ^ [z^+ir +i-^)'''kiak + hz + l/z) 

= {-l^-^k - ly.z'' 



/n-k-l\ {l-z'f 

\k-i-i) (k-ey ■ 



kzj: 



^ (z2 +mz+ 1)*^+! (z2 + 1)* 
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Proof. The definition of the gamma function as 

/•oo 

(10.12) r(s) = / t'-^e-^dt 



and the integral representation for the Hurwitz zeta function 

, 1 t"-^ e-i* dt 

(10.13) C{s,q) 



r(s) Jo 1 - e- 
are used in 



h(z) = i-l f / tfc-ie-*(-+i/-) dt - (-1)'' / ^_e-*(-+i/-) dt 
Jo 7o 1-e * 

to obtain the result. □ 

The integrals Ik{z) are now expressed in terms of the Bernoulli numbers. The 
proof is similar to the one given for Lemma |10.2[ so the details are omitted. 

Proposition 10.6. The identity 

(10.14) 

/ \ fe+l / , oo 



z2 + i; \2 ^ k^2i\ k J \z^ + l 



=1 



holds. 



According to (110. 6p . the n-th derivative of W{z) is obtained by multiplying 
(|10.10P and (110. 14p and summing over k. The coefficients V2n are then found as 

^(2n)(0) (-1)"-! 

In order to find explicit formulas for ^^'^^"^(0), (|10.4I) and (|10.14p are expanded 
in powers of z, and then the constant term in the sum is selected. Note that (|10.4p 
is of order 2;-"-'= as z — )■ 0, while (|10.14p is of order 2*^+^. So the product is of order 
Since T/F(")(z) is bounded as 2 — > 0, after summing over k all coefficients 
of z^ for i < must vanish. 

The computations to derive V2n with this approach are trivial but lengthy, and 
the resulting expression (involving multiple nested sums of binomial coefficients) is 
not particularly illuminating, so they are omitted. The vanishing of the coefficients 
of negative powers comparing it with (|6.13p yields a family of identities. 

Proposition 10.7. Let 

A{i,j,k,e,m,r) 

= / ly+J+fc (''\ fk + 2i\/2m-k-l\/k + t + m-r^j-l\ 



Then 



2m k rn rn — r—j 



*/ J>0, 
2s 



20 



M. COFFEY FT AL. 



11. Calculation of Vn by Hoppe's formula 
The function V{z) in (|4.2|) can be written as 

(11.1) V{z)^F{giz)) + \ogiz^ + l) 
with 

(11.2) F(z) = V ( - ) +log2 and 5(z) 



z2 + 1 

The expansion 

(11.3) log(z2 + 1) = y ^-^ 

n=l 

is elementary, therefore the coefficients Vn in the expansion (j4.3p are now evaluated 
from F{g(z)). 

Hoppe's formula for the derivative of compositions of functions is stated in the 
next theorem. See [12] for details. 

Theorem 11.1. Assume that all derivatives of g and F exist, then 

(11.4) f-^)"F(5(z)) = E^^^FW(g(z)), 



wh' 



ere 



(11.5) PnAgi^)) = E(-i)'-^ Q) [gi^t^' (^)" W)Y 

and P„,o(0) = /or n > 0. 

Hoppe's formula is now used to compute the n-th derivative of the function 
F{g{z)), where F is defined in (lll.2p and g{z) = z/(z^ + 1). The formula requires 

(11.6) F(^Hz)=.(^^y F{z) and (^^J [g{zW . 
These terms are computed next. 

Lemma 11.2. Let F{z) ^ ij{l/z) + log z and i/jriz) ^ {-^Y tpiz). Then, if k> I, 



r—l \ / \ 



(11.7) 

Proof. Hoppe's formula gives 



dz 



(-l)*-i(/c-l)! 



dz J \ z J ^—^ r! ' \ z J \ dz 

r=0 
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with 
(11.9) 



Pi 



Z 



k r 



^'^y^j^i-^Y H(^+fc-i)! 



yr+k 



l=Q 



e\ir - ey.ii - ly. 



^r+k -U-l 



for r > 1 and Pkfi{l/z) — 0. The last step uses the evaluation 
(11.10) 



^ i-iYr\{i + k-iy. ^ A - 1 



£1 (r -£)!(£- 1)! 



1 



Lemma 11.3. For g{z) = zj {z^ + 1) and n, j e N; 



dz 



[g{zW = n!^(-l) 
Proof. The binomial theorem gives 
(11.11) 



^ /^j + r - 1^| /^2r + ,2,-+j-n 



Differentiating n times yields the result. 

The terms in Theorem 111.11 are now written as 



k (k-1 



(11.12) F'-^Hgiz)) = (-l)'=fc!^ 



di) (z2 + l)'=+'- , (z^ + l 



I 2 , 1 \k 

H-lf-\k-iy } \ > , forfc>l, 



and 

(11.13) P„,fe(5(z)) ^z'^- 
The sum 



□ 



□ 



2n 



with F('')(g(z)) and P2n,k{g{z)) given by (|11.12p and (|11.13p . is expanded in powers 
of z. The constant term gives an expression for V2n- 
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12. An alternative approach to the valuations of Vn 

The result of Theorem 11.21 is discussed here starting from a recurrence for z„ = 
^nv2n- Using Legendre inverse relations found in Table 2.5 of [IS], the formula 
(|6.13p for V2n, namely 



(12.1) V2n = (-1) 



n+1 



+ k-l\ B2k 

n-k J 2k 



1 

n ^-^ 

k=l 

is inverted to express the Bernoulli numbers in terms of the coefficients u„. The 
authors wish to thank M. Rogers who pointed us to this inversion in [IB] . 



Lemma 12.1. // 

(12.2) £ - E 

then 



k-l\b 



2n ^-^ \ n — k J 2k 

k=l 



(12.3) fe„ = 5](-l)"-^' 



,n + k 

k=l 



ak- 



The inversion formula is used next to obtain a recurrence for a slight modification 
of the coefficients V2n- 

Theorem 12.2. Define Zn — 4nw2n — ~8ni32n(~l)- Then z„ satisfies the recur- 
rence 

r, \ n-1 / „ 

2n\ / 2n 



(12.4) -" = 2( J-E(„^J..-25.„. 



Proof. The inversion formula in Lemma ll2.1l is used with 

(12.5) a„ = 2n (^(-l)"+i«2n - ^) and 6„ = (-1)"B2. 

to obtain from Theorem 16.41 the relation 



2k- 



The result follows from here. □ 

The classical von Staudt-Clausen theorem shows that 2i?2n is a rational num- 
ber with odd denominator. The recurrence (|12.4p shows the same is valid for z„. 
Therefore 

Zn reduced modulo 2 = numerator of z„ reduced modulo 2. 

Proposition 12.3. The sequence Zn reduced modulo 2 is periodic with basic period 
{1, 1, 0}. 

Proof. The proof is by induction on n. The induction hypothesis is that the pattern 
{1, 1, 0} repeats from 1 to n — 1. 

Reduce the recurrence (|12.4p modulo 2 to obtain 

fe=l,2mod3 ^ ^ 
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This may be written as 

['~3~\ / 2 \ / 2 

k=l ^ ^ k=\ ^ 

The proof is divided in three cases according to the residue of n modulo 3. 
Case 1. Assume n — 3m. Then p2.8l) becomes 

\3m + 3fc - 2 j ~ ^ V3m + 3fc - 1 j " 

- -±( Vi. 

^ \3m + 3k-lJ 
fc=i ^ ^ 

The symmetry of the binomial coefficients shows that 

^/ 6m +1 \1 ^ / 6m + 1 \1y^/^ 6m + 1 
^\3m + 3k-l) ~ 2 ^ V3m + 3fc-l/~2 ^ Um + 3fc - 1 / ' 

k—l ^ ^ k— — ra+\ ' fc— — oo ^ 

since the terms added to form the last sum actually vanish. 
The evaluation of the sum 

may be achieved by using the WZ-technology as developed in [14:. The authors 
have used the implementation of this algorithm provided by Peter Paule at RISC. 
The algorithm shows that F(m) satisfies the recurrence 

(12.10) - 64F(m) + 65F(m + 1) - F{m + 2) = 0. 
The initial conditions F{\) = 42 and F{2) = 2730 give 

(12.11) F(m) = ^(64" - 1). 



3 



Therefore 



fc=l 



and then 

(12.13) = -i(64'" + 2) = 0mod2. 

3 

This completes the induction step in the case n = mod 3. The other two cases, 
n = 1, 2 mod 3, are treated by a similar procedure. The induction step is complete. 

□ 

Corollary 12.4. If n= 1, 2 mod 3, then V2{zn) ~ 0. 

Proof. The previous theorem shows that the numbers Z3,n+i and have odd 

numerators. □ 
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Note 12.5. The method used to obtam the values of z„ modulo 2 does not ex- 
tend directly to modulo 4 and 8. The corresponding binomial sums satisfy similar 
recurrences, but now there are boundary terms and lack of symmetry prevents the 
WZ-method to be used effectively. 
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